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MM MM TTTTTTTTTT HH HH GGGGGGGG AAAAAA TITTITITTT AAARAA NN NN 
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Mm MM TT HH HH GGGGGG AA Aa TT AA AA NN NN eeee 
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MIHSGATAND2 = Standard G tearing. ce, Tangent With 2 Arguments 
MTHSGATAND_R7 = Special GATAND routine 
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1 
1 
oating Point Arc Tangent Functions 
AN, GATAN2,GATAND,GATAND2) 

: MTHGATAN.MAR EDIT: RNH2005 


-TITLE MTHSGATAN 


l 
; T 
IDENT /2-005/ ; File 


LARRRAAAAAAAAAALAL EARLE ALAS ALALAL LAER E ARERR AAS AEE SERRA RARER SEAR SRR ERS ARSE SE SO 


COPYRIGHT (c) 1978, 1980, 1982, 1984 BY 
DIGITAL EQUIPMENT CORPORATION, MAYNARD, MASSACHUSETTS. 
RIGHTS RESERVED. 


* 

® 

® 

& 

* 

® 
FTWARE IS FURNISHED UNDER A LICENSE AND MAY BE USED AND COPIED * 
ACCORDANCE WITH THE TERMS OF SUCH LICENSE AND WITH THE ®* 
OF THE ABOVE COPYRIGHT NOTICE. THIS SOFTWARE OR ANY OTHER * 
REOF MAY NOT BE PROVIDED OR OTHERWISE MADE AVAILABLE TO ANY * 
SON. NO TITLE TO AND OWNERSHIP OF THE SOFTWARE IS HEREBY * 
® 

® 

® 

*® 

® 

es 

w 

® 

& 

® 

* 


AND 


;* TRANSFERRED. 


THE INFORMATION IN THIS SOFTWARE IS SUBJECT TO CHANGE WITHOUT NOTICE 
ND OULD NOT BE CONSTRUED AS A COMMITMENT BY DIGITAL EQUIPMENT 


OOOQOQCOooooooooooooooo 
COOQoooooooooooooooooo 


A 
CORPORATION. 


DIGITAL ASSUMES NO RESPONSIBILITY FOR THE USE OR RELIABILITY OF ITS 
SOFTWARE ON EQUIPMENT WHICH IS NOT SUPPLIED BY DIGITAL. 


* 
® 
® 
& 
® 
® 
® 
* 
® 
® 
:* OTH 
* 
* 
® 
e 
® 
® 
® 
® 
® 
® 
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; tangent value in radians of its G floating point argument. 

; MTHSGATAN2 is two argument G floating arctangent. The call is 

; standard call-by-reference. 

MTHSGATAN_R? is a special routine which is the same as MTHSGATAN 
; except a faster non-standard JSB call is used with the argument in 
; RO and no registers are saved. 


; MTHSGATAND is a function which returns the G floating point arc- 
tangent value in radians of its G floating point argument. 
; MTHSGATAND2 is two argument G floating arctangent. The call is 


BOO SOO SS OSS SS SSS OSSOS 


; standard cal “ayarererence, 
; MTHSGATAND_R S a special routine which is the same as MTHSGATAND 


; except a faster non-standard JSB call is used with the argument in 


OOoOoOooooOooooooooooooooosoooo 
SOOGOOoOSCOOCOOOOCOOOOCSOOOSOOOSOOSOOSOOOOSOOOOOOSOOOOOOOOOOSoOO 


SSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSS 


RO and no registers are saved. 
3; VERSION: 1 


; HISTORY: 
AUTHOR: 
: Steven B. Lionel, 15-Jan-79: Version 1 


4 

10 

V 

18 

14 

1? 

i 

15 

3 

4 ; 

8 ; 

4 : FACILITY: MATH LIBRARY 
1 : ABSTRACT: 

5 : MTHSGATAN is a function which returns the G floating point arc- 
5: 

8 ; 

9 ; 

40 e 

nF 

rk ; 

44; 

re 

ie 

29 

0 ; 
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: ; VERSION: 2 
: HISTORY: 
: AUTHOR: 
: Bob Hanek, 08-Jun-81: Version 2 
28 ; MODIFIED BY: 

00 
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ONO 
Sete 


oe 
MTHSGATAN ; G Floating Point Arc Tangent Functions 16-SEP-1984 01:25:15 VAX/VMS Macro V04-00 Pa 
ats ° ° aioe 94333533 ERTHRTL. SREINTHGATAN.MAR: 1 3) 


ALGORITHMIC DIFFERENCES FROM FP=11/C ROUTINE: 
1.° To avoid various flags subroutine calls have been used. 


Edit history for Version 1 of MTHSGATAN 


1-001 = Adapted from MTHSDATAN version 1-001. SBL 15-Jan-79 
1-002 = Added degree entry points. RNH 15-MAR-1981 
Edit history for Version 2 of MTHSGATAN 


-005 = Change MTHS$SAB_ATAN to MTH$S$AB_ATAN_V. RNH 29-Sep-81 

Change MTHSGATAN2D entry to MTASGATAND2 in order to conform 
with original specification. RNH 05-Oct-81 
Un-did previous edit to conform with PL/1. 
Modified small © aaa logic to avoid a microcode bug in the 
FPA. RNH 18-Dec-8 


3-00 - Use G* addressing for externals. SBL 24-August-1981 
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~SBTTL DECLARATIONS 3; Declarative Part of Module 


MTHSGATAN : 
atts BECL 


>on 


oooooooo SS SSS SS SSS SSS SS SSS SSeS Sees D— 


: INCLUDE FILES: MTHJACKET.MAR, MTHATAN.MAR 
: EXTERNAL SYMBOLS: 
. -DSABL GBL 
"EXTRN MTHSK_INVARGMAT 
EXTRN MTHSSSIGNAL 3 Signal SEVERE error 


TEXTRN MTHSSAB_ATAN_V : Gobal table used by all Arctangent 
é routines. Part of MTHATAN.MAR 


; EQUATED SYMBOLS: 


000040FC ACMASK = “M<IV, R2, R3, R4, RS, R6, R7> ; . ENTRY register mask, int 
3; ovf enabled 
: MACROS: none 
: PSECT DECLARATIONS: 
0000 -PSECT _MTHSCODE PIC,SHR,LONG, EXE ,NOWRT 


3; program section for math routines 
OWN STORAGE: none 


3 CONSTANTS: 


6_M1.0: 


00000000 0000C010 -G_FLOATING -1.0 


Sssssssssooseoo OOOOOOOCOCOOCOOSOSOOOOSOSOSO FO 
SOoOOCOCOCOOCSOOCOCCOOOCOOSOOSOOOOSOOOOOOOOSOOOSoSOoOOoO -e 


BOOCOCOCOCOCOCCOCOCOOSOOCOOOCOOCOOSOOOOOOOOOCOSOOoOo rt 
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: PARRA RAAAASAAALALA AAA ASL A SG Constants for GATAN RETAAAAAATAARAAAATAATAARAEAEAEEE 
132 ; 
1335; Each entry of the GATAN_TABLE contains the the values of XHI, GATAN_XHI_LO 
134 ; and GATA Wt HI respectively. The table is indexed by a pointer obtained 
135 ; from the ATHSSAB_ATAN_V table. The MTHSSAB_ATAN table is common to all of 
1 § 3 the arctangent Foutines and is included as part of the MTHATAN module. 
137 ; NOTE: For performance reasons it is important to have the GATAN_TABLE 
! 8 : lLongword aligned. 
140 ; 
i eALIGN LONG 
1435 GATAN_TABLE: 
144 ; Entry 0 
0000C000 FFOF3FDA 145 -QUAD “XOOOOCOOOFFOF3FDA 3 0 Oot eee ea Oe T TE TOO 
A7116ADD BF 4 18 -QUAD “XA7116ADDBF rarer 3 Mee £36286 9 a ze <4 
S66FD64F ESOCSFDA 135 Sat ew *X566F D6O4FESOCSFDA 3: 0.10506611091781236E+00 
: «Entr 
0 7F6 reg 149 Yuan a sg st 7F6 reo 3 + 1790888894637 559E+00 
DFSEBFD aay C 150 -QUAD “XDF4EBFD Erbe C : 0.1 3 5117142230E-16 
SBBSFEGF7 68455FE0 138 al .Q ao “X9BSF EGF 768455F EO 3 0.12818216111847078E+00 
: ntry 
FEC73F 15 -QUAD “XOQOOOCOOOFEC73F : 0.15621277689933777E+0 
Opes iat? shee 83 154 » QUAD “1788350 7562383 3 0: 5c0048s838027¢096-17 
66021962 DSBESFE 133 : ee *X66021962D5BE SFE : 0.15496040572616357E+00 
8000 4343 ree 139 : YOuAD *x00008000F C9D res s 8-197 9920312 0608436 +00 
A36F FC573C 158 -QUAD “X5A ob eBAr C37 C 3; 0.104 abe 59310E-16 
59EE355B AD2CSFES 122 ‘eh ae “*XS9EESSSBAD2CSFES 3 0.19278481107058049E+00 
: ntry 
FS7A3F EF 161 -QUAD “X00002000F 87A3F EF 3 0.24977041780948639E+0 
atte: A7 55 16 -QUAD “X?7 oes 99 A7D C55 3 74074643 8 rieett 360 
78054986 $h61 FEF 168 oat ae “x7 354986 461 SEF 3 Onde tess reion 4e3sae +00 
: ntry 
O0000E FB723FF3 193 -QUAD “XOOOOEQOOFB723FF3 ; =: g7eabessora3baseE-18 
196AA F CD6BC9F 1 -QUAD “X196AA37EF CO6BC9IF 3 7.27744 07258262E-16 
13D2AASS ES) FF3 i oat ary *X15D7AA3 eee ates 5 6. 6563197 96809519600 
: ntry 
9999C 800 £30) FF8 18 -QUAD “xO AL FD9OST $CBD 3 0: 12801 ra7esbcere43-16 
1321B0F 9F 3C8 170 -QUAD “X1321BOFD909F 3C8 : 0.12821747458427643E-16 
ase DEF CACE FF? 173 a -QU80 “X1A CODER CACE FF? ; 0.371753258 391603 E+00 
: ntry 
C 6763F FF 17 » QUAD “HAFSUEC9| BF OBBCBS 3 +: 9549638570591 9e946-17 
AF44EC51 8FOBBC85 174 eQUAD = “XAF44E FOBBC85 3 7.934 919294E-1 
S601FS15 br opgcts 178 ies me ¥ *X6601F61397F SSFFD 3 07463399 26837554008 900 
: ntry 
S006 DF 104004 17 » QUAD “HDEI7T 998RF FDBCAS 3 «$2827 86069161 ¢48tE-16 
171 F7DBCA 17 eQUAD = “XxDC1 6EF 7DBCA 3 2. 4916144354E-16 
s rieee aay 173 a QUAD “x96871 By e003 3 6. : $483 2068760 3€+00 
; ntry 
C0E24008 181 -QUAD *x00004000C0E24008 : 0.867 PPRSC0T OR E+00 
4 $000 § A 1 -QUAD *XC4 3 F £0983 A 3 0:36 404 468E-16 
Serhan Best \ RU RGGI: SRRSSSRSHES 
3; Entr 
00000000 89784012 185 Y QUAD *x0000000089784012 ; 0.11702833175659180E+01 
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MTHSGATAN lo b- 
yt beSLARa Declarative Part of Modul 6-SEP-1 1: MTHRTL.SR RCIM 1H ATAN.MAR; 1 
DESAE193 4FGA3CA 1 1 QUAD “XDESAE1934F4A3CA4G 52317764 4 
6CCAG04GG Aséczb08 enersthe NOE RAE atzA 6C4008 ; ‘ft: baseo9 58 599 HH 708 
O1A 01 1 ; -QUAD *x000060 O1A : 0.16449559926986694E+01 
Bee 39e8s bast A 1 1 QUAD “134 9008 es taaaascas : 0:98 ; $$ 258 17 
O3a1884 aif 1 1 -QUAD 0 Aiba : Vosesbegs 36 4 $iE 04 
0000A000 91004024 i 13 : Entry it *x0000A00091004024 ; 0 3570 0197 428955E+01 
an B4 01 194 “QUAD *XD4B oF 4 i. 18 p ; ° 
RSSGRGRD S2OB@SES SSR 135 AUAS | CHRGEARESWRDOBISTS Eo: fbouS2rBebatScteod 
RrODBCBA Ot 132 ua mrstartanyyttey i = SesaResese2 Bousee 46 
siete Beep esee 193 SQUAD  “XFBF746462E5A4016 + 0.1386 16881 417541E+01 
1 13 
1 3; Tables to be used in POLYG for couput ing GATAN: GATANTAB1 is obtained 
1 : from Hart et. al. (No. 4904). GATANTAB2 is the same as GATANTAB1 except 
1 4 ; that CO is set to 0 
ie: 
1 GATANTAB1: 
F61 1F943FD3 01 a -QUAD  *X696A2F611F943FD3 3 ¢ 0.74700604980000002E-01 
oreas 9A 43€ 38F07 1 09 ; <A TeAZOADGSESBED? ; 33 ° =: 90879628818499956 01 
¢ 71 aarp 1 10 *XBAB6DA1C71033FD 3 (4 = 0.111109) 5300320E+00 
dD BODE 4924 BFE 1 11 -QUAD “*X3DD Bor .9scer eo : C3 = 5 19882 14 198848 6E+00 
029999 2? ere 1 13 : “¥04029999999935FE 3; C2 = 0.19 39999 99893 8E+00 
5445555 SSSSBFFS 01 1 QUAD *X554A55555555BFF5 3 C1 = =.33333333333333270E+00 
00000000 9004010 1 14 QUAD * 6000000900004010 : CO = 0.10000000000000000E +01 
000000 | i GATANLEN1 = .- GATANTAB1/ A 
1 17 GATANTAB2: 
696A2F61 1F943FD 1 18 QU *X696A2F611F943FD3 3 C6 = 0.74700604980000002E-01 
974A29A9 aS 3erD7 1 19 -QUAD 4X 19 eAS9A9456 $807 : C5 = -.908796 abet oogoee Ot 
AB6DAIC 71C33FD 1 0 .QUAD DAIC71C33FDC : C4 = 0.11111091685300320E+00 
D F 49 Aare 1 1 . QUAD “K DD89 DF eaeueres : (3 = 5° 16883 14219884826E+00 
FE 1 é QUAD  *X0402999999993FE : C2 0.19 3399 3 8E+00 
5445555 5555BFF 1 QUAD *X554A55555555BFF5 3 C1 = -.33333333 270E+00 
00000000 900000 1 4 QUAD *x00000 9900000000 : CO = 0.00000000000000000E +00 
00000 GATANLEN2 = .- GATANTAB2/ 
1 G_PI: 
20185444 21FB4029 ; apn Ps *X2018544421FB4029 3 pi 
20185444 21FB4019 S gi -ovk *X2D18544421FB4019 3 pi/2 
20185444 21FBC019 | é pi ovek aap . *¥2D18544421FBC019 3 ~pi/2 
20185444 21FB4019 i 4 cnr eek y a. *X2D18544421FB4019 3 High order bits of pi/2 
5C073314 A6263(B1 $ hey RAGKD " ®¥5C073314A6263(B1 ; Low order bits of pi/2 
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; Declarative Part of Modu -SEP-1984 MTHRTL. SRE SATH HGATAN.MAR; 1 (6) 


eeerereereerereeeerereeee Constants for ATAND *teeeeerenereeeeeeereeereretet 


: ie GATAND i} My ATARSV able table is indexed by a pointer obtained 
: iron, the ATHSSAB_ATAN_V table. Theat “uth S$AB_ATAN_V table is common 
ncluded 3s part of the MTHATAN cone. 


3 ot. th °, ersten ent routines and is stu 
; : Py ae. pine: reasons it is important to have the GATAN_TAB 
rot aligned 


OWODNAUIS WOO 
Sete 


5 

& 

4 

4 

4 

4 

4 

& 

4 

4 

4 

‘ co om 
RE PSG pS CHEST, | SSE 
3 entry *x000000007F 6E SFE 128 884544372559E+0 
Snes GaSe IS tC) UCTS 
C 73FE 61° -QUAD *XQOOOCQOOFEC73FE3 ; 0.15621277689933777E +00 
SS UTA Se: STALE 
FCOD3FE 65° QUAD —*X00008000F C9D3FE ; 0.19520920515060425E+0 
i ee ca Bosc 
FS7A3F EF : ~QUAD — *X00002000F 87A3F EF ; 0.24977041780948639E +00 
cous Cote ay SHBVSateBScveSee |: TESSHESUEH SHES 
F3 + Be QUAD “X00 723FF3 ; 0.31222221255302429€ +0 
Shae SSE Bette URES BSE 
ESDISEFS 1 ~ QUAD F3DI3FFB ; 0.38988155126571655E+00 
SE ES is SHEARS Gc HASSMSUPELSEEAS 
SFFF F ~QUAD G7OSEEF : 0.49844139814376831E+00 
SHS a (STG: HEIRS 
‘<< - QUAD DF 1D40 ; 0.65223568677 E+00 
cS ISR ea eeccec Ue 
; QUAD “X COE240 ; 0.86729538440704346E+ 

iia RE as acai | GSS Meht a 
9784 $ -QUAD “xX 78401 ; 0.11702833175659180€ +01 
CPR REGS Gb SR CER: ESE 
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RATIONS ; Declarative Part of Modu MTHRTL.S SATHGATAN. MAR; 1 
D 
Bipeoea) ABiD acts 
oo ibz 


, 8 or ls 1BDG01A + 0.1644 9 
cae angorrpere if cae 
Sara 


Ww 


64495 4E+01 
7 8 Ssaehass 000 Ose 
24 ; -QUAD  *X 0004 91004024 0.25708019733428955E+01 
ree an ee ee ce ae 
Q cOuRSBe} SEARS ‘ “QUAD “HC0BB58E1 (GAABD 18 S387 so6es 3674 reset 
ECD78FEF 0B86407 ‘ - QUAD “XECDS FEF DB86407 0. 79430088028242025E+02 
10 ; iT Tables fe be used in POLYG for computing GATAND: GATANDTAB1 is obtained 
11 ni Ay outs ag ty Shs Ae coefficients given in Hart et. al. (No. 4904) by 
ig: iBo7o' is the ous as GATANDTAB1 except that "CO is set to 
173 7b - 64 eee of 180/pi 
1g GATANDTAB1: 

50669 1ECOED 1 1 -QUAD *X96450 orc Oeee) : C6 = 0.4 Say pee eT Lae: 9E+01 

AGE 1D5AC 4 FCC 8 -QUAD “*XAGEIDSACDSFCC $6 3 C5 = -.520701917 ¢ 74786E+01 
‘ 44313 6F 9 1 -QUAD mt 99E4D76F 94039 3; (4 = 8-8 36109833 060939E+01 
TeA7 bs 5EC6C040 0 QUAD “Xx 3 3 HT oo 3; C3 = -.818 Lees eeeer7ee8) 
Bie § 1 QUAD =*X/BA B8 Beier $$ ; C2 = 0.1145915 9 555564E+02 
tf sp 4 30 £33 § -QUAD X2BF Sore Hate 3 C1 = =-.190985931 Tet ef 
att A63 ASDC4 rat) -QUAD “*XC1FB1A63A5DC406C : CO = 0.57295779513082323E+02 
0000000 : GATANDLEN1 = .- GATANDTAB1/8 
§ GATANDTAB2: 

50669 1ECDES 1 -QUAD *XxX96450 precess 1 3 C6 = 0.42800293924279389E+01 
AGE1D5AC roc? 8 QUAD Retr 3 cD FCCO34 3 C5 = =-.520701917 ¢ 747 44 
SEB99E4D 76F 94039 S eQUAD “*X5E899E4D76F 94039 3; C4 = 0.63661865935060959E+01 
At 3 pecornee 0 -QUAD “Xx oes Bees eres 3 C3 = ree opty 313962579801 

8 ae ice PL 1 “ou “Xx BA7/B8 ee iecnes 3 ¢ * 0.11459155902555564E+0 
2BF066ED 193DC05 § a oe QUAD B as, X2BF 066ED193DC05 : C1 = =-.19098593171027403E+0 
FO422CE1 DIIFCO3SA . <—- * eiertaae ; CO = -.67042204869176789E+01 

00000007 5 GATANDLEN2 = .- GATAN AB2/8 
8 G_90: 
00000000 80004076 : -QUAD %*x0000000080004076 ; 90. 


40 G_M90: 
00000000 8000C076 rf 6 180: -QUAD %x000000008000C076 : 90. 
00000000 80004086 43 ~  — .QUAD *x0000000080004086 ; 180 
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RTHSGATAN » GFL 
st te ATHSGAT 
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; FUNCTIONAL DESCRIPTION: 

: GATAN = G floating point function 

; GATAN is computed using the following steps: 


1. If X > 11 then 
a. Let W = 1/X. 
b. Compute GATAN(W) = WeP(Wee2), where P is a polynomial of 


egree 6. 
set GATAN(X) = pi/2 = GATAN(W) 
/32 =< X =< 11 then 
Obtain XHI by table look 
zZ = (Xx 


“up. 
Compute = XHI)/(1 + heXHID « 
Compute GATAN(Z) = Z*P(Z2**2), where P is a polynomial of 


egree 6. 
d. Obtain GATAN(XHI) by table look-up. GATAN(XHI) will have 
two parts - the high order bits, GATAN_XHI_HI, and the low 
order bits, GATAN_XHI_LO. 
e. Compute GATAN(X) = GATAN_XHI_HI + (GATAN_XHI_LO + GATAN(Z)). 
=< X < 3/32 then 
a. Compute GATAN(X) = X + X#Q(X##2), where Q is a polynomial 
of degree 6. 
4. If X < 0 then 
a. Compute Y = GATAN(!X!) using steps 1 to 3. 
b. Set GATAN(X) = -Y, 


WWAWAWAAWAAWAANA ANIA A AAAI 
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; CALLING SEQUENCE: 
Arctangent.wg.v = MTHSGATAN(x.rg.r) 
; INPUT PARAMETERS: 


LONG = 4 ; define longword multiplier 
x = 1 * LONG ; x is an angle in radians 


IMPLICIT INPUTS: none 
OUTPUT PARAMETERS: 
VALUE: G floating arctangent angle of the argument 
IMPLICIT OUTPUTS: none 
SIDE EFFECTS: 
Signals: none 
NOTE: This procedure disables floating point underflow, enable integer 


99 ; overflow, causes no floating overflow or other arithmetic traps, an 
rts ; preserves enables across the call. 


00000004 
00000004 
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tts ATHSGATAN standard G Floating Arc Tang ~$Fb=13be 94:58:33 MTHRTL.SRCJMTHGATAN.MAR; 1 ’ (9) 

4 j--- 

d 

405 eENTRY MTHSGATAN, ACMASK ; standard call-by-reference entry 

4 $ 3; disable DV (and FU), enable IV 

4 MTHSFLAG_JACKET ; flag that this is a jacket procedure in 


6D O0000000'GF 9€ MOVAB G*MTHSSJACKET_HND, (FP) 


; set handler address to jacket 


= 

oO 

nn 

o 
SGOOCCOOCCOOOOOOoooo > 


WAW WANN AWIAIWIWIWIWAw 100 
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3 handler 
408 ; case of an error in special JSB routine 
50 04 BC SOFD 409 MOVG ax(AP), RO ; RO/R1 = arg 
6F © 10 410 BSBB  MTHSGATAN_R7 : call special GATAN routine 
04 $3 RET 3; return - result in RO 


a 7 

MTHSGATAN ; G Floating Point Arc Tangent Functions 16-SEP-1984 01:25:15 VAX/VMS Macro V04-00 Pa 11 
atts ATHSGA Ane : andard G Pigeting Arctang = $FE= 138s 94:93:33 YATHRTL . SREIRTHGATAN. MARS 1 ee (8) 

4 oak -SBTTL MTHSGATAN2 = Standard G Floating Arctangent With 2 Arguments 

p ; FUNCTIONAL DESCRIPTION: 

4 ; GATAN2 = G floating point function 

43 ; GATAN2(X,Y) is computed as following: 

D9 If Y = 0 or X/¥ > 2%#57, GATAN2(X,Y) = PI/2 * (sign X) 

D9 If Y > 0 and X/Y =< 2057, GATAN2(X,Y) = GATAN(X/Y) 

If ¥ < 0 and X/Y =< 2#57, GATAN2(X,Y) = PI * (sign X) + GATAN(X/Y) 


; CALLING SEQUENCE: 
Arctangent2.wg.v = MTHSGATAN2(x.rg.r, y.rg.r) 


; INPUT PARAMETERS: 
00000004 x = 1 * LONG ; x is the first argument 
00000008 y = 2 * LONG : y is the second argument 


; SIDE EFFECTS: See description of MTHSGATAN 


40FC ENTRY MTHSGATAN2, ACMASK 3 standard call-by-reference entry 
3; disable DV (and FU), enable IV 
MTHSFLAG_JACKET ; flag that this is a jacket procedure in 


6D Q0000000'GF 9 MOVAB G*MTHSSJACKET_HND, (FP) 


3D 
3D 
3D 
3p 
03D 
03D 
03D 
03D 
O3E ; set handler address to jacket 
O3E 3: handler 
O3E 444 3; case of an error in special JSB routine 
50 04 BC SOFD bae$ 445 MOVG ax(AP), RO ; RO/R1 = argl 
52 08 BC SOFD tr rr ; MOVG @y(AP), R2 3; R2/R3 = arg2 
0 4 33 : Test if Y = 0 or X/Y > 2#*57 
33 13 O3EC 450° BEQL. _—INF ; branch to INF if Y = 0 
54 23 OOF 8F AB 4 451 BICW pa ngOOr RO, R4 3: R4 = exponent (xX) 
55 2 OOF 8F F $36 BICW #*x800F, R2, RS : RS = exponent(Y) 
4 A FA 4 SUBW R5, R4 : R4& = exponent(X) = exponent (Y) 
O3A0 8F 4 61 FD 454 CMPW R4, #58*16 3 compare R4 with 58 
14 40¢ 439 BGTR INF : if K/Y > 2857, branch to INF 
4 43 : Test if Y>Oory<0 
52 BS 0404 138 ’ TStw oR 3 test the sign of Y 
16 14 $ 460 BGTR A2PLUS ; branch to AZPLUS if Y > 0 
50 83 4 461 TSTwW 3; test the sign of X 
09. ri "1 $66 BGEQ A1PLUS : branch to ATPLUS if X >= 0 
49¢ 464 : Y < 0 and X < 0 and X/Y =< 20857 
40C 465; 
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tts ATHSGATAN2 : tandard G Pigeting Arctang ~$EE=138e 9:53:37 MTHRTL.SRCJMTHGATAN.MAR; 1 ni (8) 
35 10 04 § BSBB MTHSGATAN_R7D 3; RO/R1 = GATAN(X/Y) 
50 FDBS CF 42FD 04 SUBG2 Pl, RO 3 RO/R1 = -PI + GATAN(X/Y) 
04 04 3 RET 3 return 


; return 


P 
$ 
4 
4 
4 
4 re H 
" $f9 ; Y < 0 and X > O and X/Y =< 2**57 
4 ore A1PLUS: 
2c 10 04 47 BSBB MTHSGATAN_R7D 3 os yd = GATAN(X/Y) 
50 FDAC CF 4OFD 04 474 ADDG2 PI, RO 3; RO/R1 = PI + GATAN(X/Y) 
04 r $2? 3 otek 
4 477 ; Y > 0 and X/Y =< 2257 
4 478 ; 
4 479 A2PLUS: 
23 10 04 480 BSBB MTHSGATAN_R7D 3; RO/R1 = GATAN(X/Y) 
04 re rtf RET 3; return 
4 cas : Y = 0 or X/Y > 2e*57 
Be 484 ; 
4 489 INF : 
50 B65 04 486 TSTW RO ; test the sign of X 
08 14 3 487 BGTR 1$ : branch if xX > 0 
0¢c 13 04 488 BEQL 2s : branch if xX = 0 
50 FDAD CF 7D 04 489 MOVO G_MPI_OVER_2, RO ; RO/R1 = GATAN(X/Y) = =PI/2 
04 ‘ $29 RET 3 return 
50 FDOF CF 7D 4 G_PI_OVER_2, RO ; RO/R1 = GATAN(X/Y) = PI/2 
4 


se 
496 ; Here if both X = 0 and Y = 0. Signal INVALID ARG TO MATH LIBRARY 


OSOooo 


ere 


50 01 OF 79 2$: ASHQ #15, #1, RO RO/R1 = reserved operand, copied 


to CHFS_MCH_S SAVRO/ 1 so aanee ers 


REER 
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: can change Tf t ntinu 
7E 00" 8F oA 7 MOVZBL ple INVARGMAT, -(SP) ; code lig AR Hd Kremat to TO mATH L TBRARY 
00000000'GF 01 FB 04 CALLS , G*ATHS$$SIGNAL : Signal SEVERE error 

04 044 ; return if a handler says SS$_CONTINUE 


sey 
2-00 


50 52 46FD 
0 53FD 
0 18 
0089 


56 50 3FD8 4 ‘ 
56 006d ef B 
“2 


a. oe. eae 


56 FFFFFFOO 8F CA 
53 00000000 ' GF oF 
53 00000000' GF 


56. 6346 0 
56 FB87 CF46 = 7E 


“Ea 


5054 


% Bun 
FCB4 CF 55D 
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Special GATAN routine -SEP=1 MTHRTL.S 
~SBTTL MTHSGATAN_R7 = Special GATAN routine 
Special GATAN = used by the standard routine, and directly. 


CALLING seouentey: 
bove anything needed in RO:R7 


JSB MTHSGATA 
return with ae . RO/R1 
Note: This routine is written to avoid causing any eateger overflows, floating 


Page 1 
° ( 


IO 


; input in RO/R1 


++ hahaa or floating underflows or divide by 0 conditions, whether enabled or 
not. 
REGISTERS USED: 
RO/R1 = Floating argument then result 
RO:R5 = POLY 
R6/R7 = Y Guring POLYG 
MTHSGATAN_R7D: ; for local use only! 
dIv R2, RO 
MTHSGATAN_R?:: 3 Special GATAN routine 
TSTG RO ; R6 = X = argument 
BGEQ POS_ARG 
BRW NEG_ARG ; Branch to negative argument logic 


Argument is positive 


Bos_ARG: suays #X3F0B, RO, R6 Argument is less than 3/32 


= DOONAN WO OD NAULS WN O ODNAULS WN OVOnNG- NSW" O00 
ee 


branch to small arqumees, fogic 
CMPW #°X006D R6 : Argument is greater that 1 
BLSS LARGE _ARG ; branch to large argument tate 
: Logic for positive medium sized arguments. Get pointer into GATAN_TABLE. 
, ROTL #1, R6, R6 ; R6 = index into MTHSSAB_ATAN table 
ICL #-256 : zero high order bits of index 
MOVAL G*MTHSSAB_ATAN_V, R3_ ; RS = 4 ress of RTL vector entr 
ADDL Ht “ATAN_V, R3 ; R3 = address of MTHS$$AB ay table 
MOVB (R3)CR R6 ; R6 = offset into GATAN_TAB 
MOVAQ GATAN_ ROASLE R6], R6 : R6 = pointer to XHI 
: Compute z 

ova (R6)+ 3; R2 = XHI 

é MULG Re. R m6. > RG = XeXHI 
ADDG : RG = 1 + X*XHI 

4 SUBG R2, RO ; RO = X = XHI 

5 DIVG R4, RO 3; RO = Z = (X = XHID/(1 + X*XHI) 

5 : Evaluate Z*P(Z**2) 

i : mova (SP) : P sh Z onto the stack 

? MULG2 RO, R Ze%2 
POLYG RO, #GATANLENI-1, catantia| 

2 3 RO = P(Z**2) 


-- = 


MTH 
Tab 


C 
MTHSGATAN ; G Floating Point Arc Tangent pelet ions 16-SEP-1984 01:25:15 VAX/VMS Macro Vv04-00 P 14 
B 008 ATHSGATAN_R? = Special GATAN routine ~SeemT8bS 91:58:57 PANANAS Bacte MOK Oe man: 2° 18, 
gf 44FD A MULG (SP)+, Re 3 RO = GATAN(Z) = Z*P(Z**#2) 
rea A ADDG (R6)+, R ; RQ = GATAN_XHI_LO + GATAN(Z) 
66 40FD A ADOG (R6), RO ; RO = GATANTX) = GATAN_XHI_HI + 
8 ; (GATAN_XHI_LO + GATAN(Z)) 
8 RSB ; Return 
SMALL: BRW SMALL ARG ; Dummy Label used to avoid adding 


s an extra insrtuction in the 
; medium argument logic 


; Large positive argument logic. 


LARGE_ARG: 


56 FB46 CF 50 47FD 04B4 DIVG3 RO, G_M1.0, R6 : RO = “WwW = -1/KX 

50 56 56 45FD 048 MULG RS, RS, R * RO = Wee2 
FC91 CF «606 «50 S5FD 4c POLYG RO, #GATANLEN1=1, GATANTAB1 

; RO = P(Wee2) 
50 56 44FD MULG2 R6, RO > RO = GATAN(W) = <WeP(Wee2) 
50 FD18 CF 40FD ADDG2 G_PI_OVER_2_L0, RO : 
50 = FDOA CF 40FD ADDG2 G-PI“OVER-2-HI. RO + RO = GATAN(X) = PI/2 = GATAN(W) 
05 $B : Return 


; Logic for negative arguments 


NEG_ARG: 


SUBW3 #*XBFD8, RO, R6 Argument is less than 3/32 


RMMMMMMMMMoOOOCOTVTVT VT VTVITINAOoO 


ot tt at at = I =I DOOD ODOOO0 000 OOO OOOO 0900 09 09 09 09 09 09 09 09 SI NINN NINN NINO AAA AAO 
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4 
: : 
4CB 
4D1 
r 7 
ocbe 
O¢08 
4 
4 
a3 
A3 04 : 
19 04 BLSS SMALL_ARG ; branch to small argument logic 
56 006D ef B1 04 CMPW #*x008D, R6 ; Argument is greater than 11, 
4 19 3 3 BLSS N_LARGE_ARG ; branch to large argument logic 
3 4 : Logic for negative medium sized arguments. Get index into GATAN_TABLE. 
56 56 FF 8F 9C 04E7 : ROTL #1, R6, R6 : R6 = index into MTHSSAB_ATAN table 
56 FFFFFFOO 8F CA 04 g BICL #-256, R6 ; clear high order (unused) bits of ind 
53 000000'G oF 4 MOVAL G*MTHSSAB_ATAN_V, RF ; R3 = address of RTL vector ye 
53 0°G GFA ADDL G*MTHSSAB_ATAN_V, R ; R3 = address of MTHSSAB_ATAN table 
6 6346 0 0501 MOVB (R35) CR6) Re : R6 = offset into GATAN_TABLE 
56 FAFE CF46 7E 3 5 MOVAQ GATAN_TABLECR6], R6 : R6 = pointer to XHI 
3 ; Compute 2 
52 : 7D 05 y MOVQ. = (R6) +, R2 ; R2 = XHI 
54 45FD 3 MULG Re, RO, R4 > RG = X*XH] 
54 98 4 43FD 1 SUBG R4, #1, R4 3; RG = 1 = X*XHI = 1 + X*(-XHI) 
2 40FD 051 ADDG R2, R 3; RO = X + XHI = X = (=XHI) 
50 4 46FD 2) DIVG R4, R ; RO=2 
5 : Evaluate Z*P(Z*#2) 
43 7D 2 mova RO, -(SP) ; Push Z onto the stack 
3 44FD : MULG2 Re. RO > RO = Ze"2 
FC2A CF 606 55FD POLYG RO. #GATANLEN1-1, GATANTABI 
52E 3; RO = P(Z*#2) 


Dd 
MTHSGATAN G Floating Point Arc Tangent lekettons 16-SEP-1984 01:25:15 VAX/VMS Macro V04-00 P 15 Tt 
Orbos WTRSGATANAY  Spechol GATAN tout tge eet aes MTMRTL«SREIMTHGATAN.MAR:1 O° (9) 2=( 
8 HH 44GFD § 0 MULG (SP)+, Re ; = GATAN(Z) = zeP(zeee) 
ria 1 SUBG (R6)+, FY = GATAN_XHI_LO + GATAN(Z) 
0 66 42FD 6 § (R6), RO 3 . = SATAN? = GATAN_XHI_HI + 
A 6 $ (GATAN_XHI_LO + GATAR(Z)) 
05 é ‘ ? : ; RSB ; Return 
: ; ° § ; Logic for large negative arguments 
538 608 * 
538 6 : N_LARGE_ARG: 
56 _ 56 36 rh 4a 348 6 ? pive Re’ G gM » R6 ; Re s —— 
Fcoa tr 0G 8G $f 4 é ; POLYG RO. SOK TRNLENT= ‘ GATANTABI 7 pieesaie 
MH s ® 
50 56 44FD 2 634 MULG R6, RO ; RO = GATAN(W) = WeP(Wee2) 
50 FC91 CF 4geD 552635 SUBG G_P1-OVER2_L0, ; 
50 «FC83 CF 42FD 0558 «= 6 6 SUBG2 G_PI-OVER~2-HI, ; RO = GATAN(X) = GATAN(W) = PI/2 
05 OSSE 6 $B ; Return 
aeF 836 ; 
O55F 640 ; Small argument logic. 
beer bee | 
035 $8 SMALL_ARG: 
56 50 7D OSSF 644 mova 3 R6 = prqument = Xx 
50 $008 F AA 268 645 BICW i: RO 3; RO = 
50 —70 8F BI 036 eae CMPW #*x : Fonpere "2%=26 to IX! 
19 Q056C 64 BLSS ; Needs polynomial evaluation 
50 56 7D $F os mMOva Ae. RO ; Return with answer equal to 
05 2 $63 RSB § argument 
50 50 44FD 258 651 1$ MULG2 RO, RO ; RO = Xe22 
FCI3 CF §606)Ss 50 SSFD 376 63¢ POLYG RO, #GATANLEN2-1, GATANTABZ esc 
50 56 44FD 057D 654 MULG2 R6, RO > RO = X#0(Xee2) 
50 56 40FD 44 655 ADDG R6, RO 3; RO = GATAN(X) = X + X#Q(X##2) 
05 0585 $26 RSB ; Return 
586 665 
nn — 


IMPLICIT INPUTS: none 
OUTPUT PARAMETERS: 
VALUE: G floating arctangent angle of the argument 
IMPLICIT OUTPUTS: none 
SIDE EFFECTS: 
Signals: none 
NOTE: This procedure disables floating point underflow, enable waheger 


; overflow, causes no floating overflow or other arithmetic traps, an 
; preserves enables across the call. 
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bos ATHSGATAND : Standard G Pigeting Arc Tan atts 94:53:33 LMTNRTL. SREINTHGATAN.MAR: 1 = 18) 
686 «SBTTL MTHSGATAND = Standard G Floating Arc Tangent 
# p++ 
? 007 : FUNCTIONAL DESCRIPTION: 
2 ? 685 ; GATAND = G floating point function 
: 6 ge ; GATAND is computed using the following steps: 
5 4 6 : > 61, If X > 11 then 
2 ? ° 3 a. Let W = 1/X, 
3 b. Compute GATAND(W) = WeP(W**2), where P is a polynomial of 
5 ? 672 ; degree 6. 
5 673 ; c¢. Set GATAND(X) = pi/2 = GATAND(W) 
586 674 ; 2. If 3/32 =< X =< 11 then 
586 «675 ; a. Obtain XHI by table Look-up. 
586 of8 ; b. Compute Z = (X = XHI)/(1_+ X*XHI). 
0386 of ; Se Compute GATAND(2) = Z2*#P(Z7**#2), where P is a polynomial of 
: egree 6. 
586 ere ; d. Obtain GATAND(XHI) by table look-up. GATAND(XHI) will have 
586 ©6680 ; two parts - the high order bits, GATAND_XHI_HI, and the low 
586 «= 681 —; order bits, GATAND_XHI_LO. 
386 6 ¢ ; e. Compute GATAND(X) = GATAND_XHI_HI + (GATAND_XHI_LO + GATAND(Z)). 
82 6 6835; 3. If 0 =< X < 3/32 then 
280 684 ; a. Compute GATAND(X) = X + X#Q(X*#2), where Q is a polynomial 
0586 685 ; of degree 6. 
be 6 O36 3 4. If X < 0 then 
586 oe 3 a. Compute Y = GATAND(!X!) using steps 1 to 3. 
0586 688 ; b. Set GATAND(X) = -Y. 
0886 690 
0586 691 ; CALLING SEQUENCE: 
086 693 ; Arctangent.wg.v = MTHSGATAND(x.rg.r) 
05 6 699 > INPUT PARAMETERS: 
00000004 6 69 LONG = 4 : define longword multiplier 
00000004 ? 843 x = 1 * LONG ; x is an angle in degrees 
6 
6 
6 
6 
6 
6 
é 
6 
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loat 
ATAN 
40FC E 
5 
i 
6D O0000000'GF 9€ 
3 F 
58F 
S8F 
50 04 BC SOFD O58F 
6F 10 0594 
04 B228 
059 


.ENTRY MTHSGATAND, ACMASK 
MTHSFLAG_JACKET 


MOVAB G*MTHSSJACKET_HND, (FP) 
MOVG  ax(AP), RO 
BSB MTHSGATAND_R7 


£ 
int Arc Tangent FunctTons 16-sEp-1984 01:25:18 
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MTHRTL. 


standard call-by-reference entry 
disable DV (and FU), enable IV 
flag that this is a jacket procedure in 


set handler address to jacket 
handler 


case of an error in special JSB routine 
RO/R1 = arg 

call special GATAND routine 

return - result in RO 


(10) 
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G 
stagger aPocaars Toppan furettone SEEECISES 1:35:35 YAMARE Meee AOtGHD uns Pome Hf] |B 
= -SBTTL MTHSGATAND2 = Standard G Floating Arctangent With 2 Arguments 
FUNCTIONAL DESCRIPTION: 
GATAND2 = G floating point function 


GATAND2(X,Y) is computed as following: 


o— 
mo 


If Y = 0 or X/Y¥ > 2%#57, GATAND2(X,Y) = 90 * koto x) 
If Y > 0 and X/Y =< 2e*57, GATANDG (X.,¥) = GATAND(X/Y) 
If ¥ < 0 and X/Y¥ =< 2257, GATAND2(X,Y) = 180 * (sign X) + GATAND(X/Y) 


CALLING SEQUENCE: 
Arctangent2.wg.v = MTHSGATAND2(x.rg.r, y.rg.r) 
INPUT PARAMETERS: 


x = 1 * LONG 3; x is the first argument 
y = 2 * LONG ; y is the second argument 


3; SIDE EFFECTS: See description of MTHSGATAND 


00000008 


40FC -ENTRY MTHSGATAND2, ACMASK 3 
MTHSFLAG_JACKET : 


MOVAB G*“MTHSSJACKET_HND, (FP) 


ard call-by-reference entry 
le DV (and FU), enable IV 
that this is a jacket procedure in 


NNN SII 
PUD ASI BF BS BS BS SB AAA AMAA AID 


NOAUALS WN  O OO NAUE WN" OOONAUS UT" OO 


6D  Q0000000'GF 9E 
3; set pentior address to jacket 


PPP PPS FP OOOO OOOOOOOOOOOOOOOVOOOOOOOOOOOOO Pro 


QOOOOoWWW> 


3; handler 
3: case of an error in special JSB routine 
50 04 BC SOFD MOVG ax(AP), RO : RO/R1 = arg 
52 08 BC S5OFD MOVG @y(AP), R2 ; R2/R3 = arg2 
; Test if Y = 0 or X/¥ > 2857 
‘ BEQL_ INF BEG ranch to INF_DEG if Y = 0 
6 29 oor BICW #*xXBOOF, RO, R4 4 = exponenttXx) 


BICW3 #*XB00F. R2. RS 
SUBW RS, -R4 

R4, #58*16 
BGTR _INF_DEG 


: Test if Y>Oory<0 


4 = exponent(X) = exponent (Y) 
epore R4 with 58 
f X/Y > 2**57, branch to INF_DEG 


O3A0 8F 


— WWD COL 


b 
R 
R = exponent (Y) 
C 
i 


SON 
NNN NIA AAA APAAAAOUIM 
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wooon CFunnw 
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TSTW R 3; test the sign of Y 
5C4 BGTR A2PLUSD : branch to AZPLUSD if Y > 0 
266 TSTW 3; test the sign of X 
ef BGEQ A1PLUSD : branch to ATPLUSD if X >= 0 
sth : Y <0 and X < 0 and X/Y =< 2#*57 
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MTHSGATAN 
2-008" 
35 
50 = FDEF CF 42 
2c 
50 = FDE6 CF 
23 
50 
08 
oc 
50 =FDCF CF 
50 FDC CF 
50 01 OF 
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SSEe=158 11:98:23 


oating Arctan 


BSBB aa a R7D p pore = GATAND(X/Y) 
SUBG2 : = -180 + GATAND(X/Y) 
RET 3; return 

: Y <0 and X > 0 and X/Y =< 2*#57 

Aipwuso: 
BSBB MTH TEGATAND. R7D : ad = GATAND(X/Y) 
ADDG2 6.18 RO $ 1 = 180 + GATAND(X/Y) 
RET 3; return 


a 
ad 
ae 


3; Y > 0 and X/Y¥ =< 257 

A2PLUSD: 
BSBB MTHSGATAND_R7D 3; RO/R1 = GATAND(X/Y) 
RET 3 return 

:; Y = 0 or X/¥ > 2**57 

INF _DEG: 
TSTW RO 3; test the sign sf x 
BGTR 1$ 3 branch if X > 
BEQL 2s 3 branch if xX = 0 
mova G_M90, RO 3; RO/R1 = GATAND(X/Y) = -90 
RET : return 
MOVva 6.90, RO ; RO/R1 = GATAND(X/Y) = 90 
RET 3; return 

Here if both X = 0 and Y = 0. Signal INVALID ARG TO MATH LIBRARY 

ASHQ #15, #1, RO 


FAP Se te ee 
of + 


RO/R1 = reserved operand co180ed 
to CHFS$ NCHS SAVRO/ 1 so handlers 
can change want to continu 
— tore Twat ID 
sognes SEVERE e 
return if a handler says SS$_CONTINUE 


oie 8 dda + INVARGMAT, =(SP) 


, G*ATHS$SIGNAL 


Page 19 


(11) 


De knGaat TO MATH LIBRARY 
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on = Special gat AND Petia -SEP-1 MTHRTL SRE SATH HGATAN.MAR; 1 
«SBTTL MTHSGATAND_R7 = Special GATAND routine 
: Special GATAND = used by the standard routine, and directly. 
, > CALLING SEQUENCES: 
: pove anyth needed in RO:R7 
§ 3 ; input in RO/R1 
; 138 MTHSGATAND_R7 
8 3 return with result in RO/R1 
3; Note: This routine is written to avoid causing any integer overflows, floating 
0; = ica or floating underflows or divide by 0 conditions, whether enabled or 
3 not. 
: REGISTERS USED: 
¢ ; apve /R1 = Flossing argument then result 
; ; Porn? = -Y y ee POLYG 
40 wmenns ). R7D: : for local use only! 
50 52 46FD 41 Vee R2, RO 
rX: RTHSGATAND R 3 Special GATAND routine 
50 S3FD 4 TSTG RO ; R6 = X = argument 
18 44 + oy POS_ARGD 
00 31 o2 BRW NEG_ARGD ; Branch to negative argument logic 
oF : Argument is positive 
rk bos _ARGD: 
56 50 3FD8 BF A3 50 SUBW3 = #*X3FD8, RO, R6 :; Argument is less than 3/32 
A 19 51 BLSS MALLD ; branch to small argument logic 
56 006D % 26 CMPW #*x006D : Argument is greater that 11, 
6 19 3 BLSS —_LARGE ORaGD. : branch to Large argument logic 
22 3 ; Logic for positive medium sized arguments. Get pointer into GATAND_TABLE. 
6 8 9C 8 ' ROTL #1, R6, R6 : R6 = index into MTH$SAB_ATAN table 
sa creer : cA 28 BICL : zero high order bits of index 
2 43 5 MOVAL GaRTHESAB. aren -Ve a 3 ns = address of RTL vector entr 
¢ 8 60 ADDL G*MTHSSAB ; R3 = address of MTHSS$AB_ATAN table 
56 o3e 6 61 MOVB (R3)CR oie ; R6 = offset into GATAND_TABLE 
56 FBB1 CF46 7E $8 MOVAQ  GATAND_ ae R6 =: R6 = pointer to XHI 
és | Compute 2 
2 7D 66 ’ MOV —- (R6)+, R2 : R2 = XHI 
54 45FD 6 MULG Re, nO, R4 : RG = X*XHI 
4 40FD rt: ADODG #1, R4 : RG = 1 + X*XHI 
4 42FD : SUBG 2, 4 ; RO = X = XHI 
46FD Y DIVG R4, R : RO = Z = (KX = KHID/(1 + X*XHI) 
he : Evaluate ZeP(zZ*#2) 
ze 7D 74° mova = RO, _-=( SP) : Push Z onto the stack 
be 44FD 75 MULG2 a : = [22 
FCDE CF 55FD 76 POLYG RO, #GATANDLENI-1, GATANDTAB1 


i  Y f 


oe 


009F 
56 = F988 gf 0 
50 6 8 

FCBB CF 8606 
50. 56 
50 FD22 CF 
56 50 BFDS BF 
56 006d HH 
4 


5 FFFFFE 
5 000 * GF 
5 * GF 
6 6346 
56 FB2E CF46 
54 


FCSA CF 


; Large positive argument logic. 


LARGE _ARGD: 
DIVG 


O00 00000009 00 C0 CO COCO OD CO II 


FV - OOONAUE US" OOO 


; RO, G_M1.0, R6 

MULG RG, Rb, R 

99 POLYG RO. #GATANDLEN1-1, GATAN 
38 MULG2 R6, RO 

98 ADDG2 G6_90, RO 

eo RSB 


901 ; 
O06 ; Logic for negative arguments 


905 NEG_ARGD: 
90 


SUBW3 - #*XBFD8, RO, R6 
90 BLSS § SMALL_ARGD 
908 CMPW 3836s #*X008D, R6 
BLSS § N_LARGE_ARGD 


Logic 


#-1, R6, Rb 
G*MTHSSAB_ATAN_V, R3 
G*MTHSSAB_ATAN_V, R 
MOVE (R3)CR6],°R 
GATAND_TABLECR6], R6 


10 
i 
13 
14 
13 
9 
18 MOVAQ 
19 ; 
9 ; Compute Z 
‘ ; VQ  (R6)+, R2 
MULG3 R2, RO, R4 
4 SUBG3 R4, #1, R4 
5 ADDG R2, R 
¥ DIVG2 R4,R 
3 ; Evaluate Z*P(Z**2) 
? , mova RO, (SP) 
MULG2 
§ POLYG R 


WODOOOOOOOOOOODODODOOOOOOOOOO 


J 

ting Point Arc Tangent Func lens 16-SEP-1984 01: 

"A? = Special GATAND routine pease 1:38 
MULGZ —(SP3+, RO 
ADDG enes R 
ADDG (R6), RO 
RSB 

SMALLD: BRW SMALL_ARGD 


3; R6 = -W = 
: RO = 
bTAB1 

; RO = 


for negative medium sized arguments. 


° 4 R = 
; #GATANDLEN1-1, GATANDTAB1 
; RO = 


; Push 


:15  VAX/VMS Macro v04-00 Page 21 mT} 
} MTHRTL.SRCIMTHGATAN.MAR:1 |” (9d) 2=( 
RO = P(Z##2) 
RO = GATAND(Z) = 2#Q(Z*#2) 
RO = GATAND_XHI_LO + GATAND(Z) 

= GATANDTX) = GATAND XHI HI ¢ 

(GATAND_XHI_LO + GATAND(Z)) 

Return 


Dummy Label used to avoid adding 
an extra insrtuction in the 
medium argument logic 


-1/x 


Wee2 


P(Wee2) 

RO = -GATAND(Z) = <Z*P(We*2) 

RO = GATAND(X) = 90 = GATAND(Z) 
Return 


Argument is less than 3/32 
branch to small argument logic 

Argument is greater than 11, 
branch to Large argument logic 


Get index into GATAND_TABLE. 
R6 = index into MTH$$AB_ATAN table 


clear high order (unused) bits of ind 
R3 = address of RTL vector entr 

R3 = address of MTHSSAB_ATAN table 

R6 = offset into GATAND_TABLE 

R6 = pointer to XHI 

R2 = XHI 

R4 = X*XHI 

R4 = 1 = X*XH] = 1 + X*(-XHI) 

Re = + XHI = X = (=XHI) 


Z onto the stack 
an2 


P(Z**2) 


x 
mT TAN ; G Floating Point Arc Tangent Func = 16-SEP-1984 01:25:15 VAX/VMS Macro Vv04-00 Page 
BOOS ATHSGATAND. A? = Special GATAND rockion pant ioets 94:33:33 MTHRTL.SRCJMTHGATAN.MAR; 1 ° (93) 
64FD £6 4 MULG2 (SP)+, R z RO = GATAND(Z) = Z2*P(Z*#2) 
agtD A SUBG (Ro), R ; RO = GATAND_XHI_LO + GATAND(Z) 
42FD 06E $ SUBG2 (R6), RO : RO = GATANDTX) = GATAND_XHI_HI + 
F : (GATAND_XHI_LO + GATANB(Z)) 
05 f : RSB + Return 
F ; ; Logic for large negative arguments 
F Bs 
F 43 N_LARGE_ARGD: 
56 F907 CF 50 47FD 06F 44 DIVG3 RO, G m1 0, R6 : RS = WwW 1/!x! 
50 56 6 $3rD FA 945 MULG RS Rb, R : RO = Wee2 
FC3A CF S5FD (O6FF 46 POLYG RO. #GATANDLEN1-1, GATANDTABT | nena 
50 56 44FD 708 348 MUL G3 R6, RO : RS = GATAND(W) = WeP(Wee2) 
50 FCA1 CF 42bD 70A 94 SUBG2 6_90, RO : RO = GATAND(X) = GATAND(W) = 90 
0 710 950 RSB 3; Return 
ay 983 
711 338 > Small argument logic. 
Hi ge 
711 996 SMALL _ARGD : 
56 30 5088 ny 32 pg RO. R6 : R6 = argument = X 
50 8000 $P AA a 336 BICW #718000, RO : RO = IX! 
50 E70 8F 261 71C 960 CMP #*x3E70, RO 3; Compare 2*=-26 to iX! 
19 0721 961 BLSS 1$ 3; Needs polynomial evaluation 
50 56 FC80 cr 45fD , 3 36¢ MULGS G,P1_0v_180_M_64, R6, RO; RO = X*(pi/180 = 64) 
50 0 44FD 07 964 1$: MULG2 RO, RO + RO = Xee#2 
FC41 CF 306 55FD , 6 369 POLYG RO. #GATANDLEN2-1, GATANDTAB2 step 
° s 
50 56 44FD 0737 367 MULG2 6, RO > RO = X*Q(X##2) 
56 0060 BF Ad 07 968 2S: ADDW #*X60, R6 > RE = Xe2ee6 
50 56 40FD 0740 969 ADDG2 R6, RO + RO = GATAND(X) = X#2e26 + X#Q(X##2) 
744 970 38: RSB 3; Return 
745 971 
745 76 
745 97 END 


MTHSGATAN 
Symbol table 


ee 


=e a 


INF _DEG 
LARGE _ARG 
LARGE —ARGD 
LONG 


MTHSSAB_ATAN_V 
MTHSSJACKET_AND 
MTHS$SIGNAL 


m 
MTHSGATAN_R7 
MTHSGATAN-R7D 
nTHSK INVARGMAT 


L 
; G Floating Point Arc Tangent ledettone = 


OOOOQOooooooooo 
= 
—=§ 90 |"WhWooounfui 
ODOC VNNLSOONG—-oO— 


Bowne 


FN" WDWOOOCWWDDOOCOCWOBVNVNODWOVNNOOMUWIMN 


SOoOoOOSOSCSOSOOOOOOOOOOOSOOOOOSOoOSoOoOO 


SOSoOOooooooooooSo 


PoP ol eee 
NWONMPMMooeo 


SOoooooo 


So 
oO 
So 
™ 


geeaeaeee 
eeeeeee 
Reeeeeee 


1 
1 
1 
1 


ooo 
—— 


OSGOOGOCOGOOOOOCGOOOOCOSSSSSOSS GCOCOCOCOOOOOOOOOoOOoOO 
th I IOS III MIOVMO 


SEP=1 
SEP=1 


3 


4 
4 


? 


1: 
1: 


33:3} 


AX/VM 
MTHRT 


S 
L. 


Mac 
SRC 


o V04-00 


r 
IMT 
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Psect synopsis ° ° g-8 sets 94:53:37 MTHRTL.SRCJMTHGATAN.MAR; 1 . (93) 


tere masa ee me moenan$ 


! Psect synopsis ! 


emma momranaaecaena $ 


PSECT name Allocation PSECT No. Attributes 
OS . 00080000 ( 9-3 0 ¢ 0.) NOPIC USR CON ABS _ LCL NOSHR NOEXE NORD NOWRT NOVEC BYTE 
_MTHSCODE 0000745 ( 1861.) O01 ¢ 1.) PIC USR CON REL LCL SHR EXE RD NOWRT NOVEC LONG 


$e eer meme wm eam aman e mer aaes $ 


! Performance indicators H 


eestor e nn eeaeeeaeen eeeee 


Phase Page faults CPU Time Elapsed Time 
Initialization 34 0:00:00.0 0:00:01.17 
Command processing 114 :00: 23 :00: 3-80 
Pass 120 : 3: 54 :00:07.51 
Symbol table sort :00:00.03 0:00: 0-33 
Pass 2 18 0:00:02.14 + Bh. 8 
Symbol table output 8: : +06 :00:00.0 
Psect synopsis output :00:00. 0: B88 8 
Cross-reference output 0:00: 2: 0 0:00: 8:8 
Assembler run totals 46 0:00:05.51 0:00:20.98 


ire working ssf Limit was 1050 pages. 

16379 bytes (32 pages) of virtual memory were used to buffer the intermediate code. 

There were 10 pages of symbol table space allocated to hold 54 non-Local ns 8 local symbols. 
1033 source Lines were read in Pass 1, producing 24 object records in Pass 2. 

1 page of virtual memory was used to define 1 macro. 


were ese mse nec eceeceee= ee cea & 


Macro Library name Macros defined 
~$255SDUA28: (SYSLIBISTARLET.MLB; 2 egy Ky 
0 GETS were required to define 0 macros. 

There were no errors, warnings or information messages. 


MACRO/EWABLE =SUPPRESSION/DI SABLE=(GLOBAL , TRACEBACK) /LIS=LIS$:MTHGATAN/OBJ=OBJ$:MTHGATAN MSRC$:MTHJACKET/UPDATE=(ENHS :MTHJACKET) +MSRC 


NT CORPORATION 
D PROPRIETARY 


AH-BT13A-SE 
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